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Abstract Type la supernovae (SNe la) play a key role in measuring cosmological param- 
eters, in which the Phillips relation is adopted. However, the origin of the relation is still 
unclear Several parameters are suggested, e.g. the relative content of carbon to oxygen 
F^ ■ (C/O) and the central density of the white dwarf (WD) at ignition. These parameters are 

Mh' mainly determined by the WD's initial mass and its cooling time, respectively. Using the 

Q . progenitor model developed by Meng & Yang, we present the distributions of the initial 

}-{ ' WD mass and the cooling time. We do not find any correlation between these parameters. 

C/3 , However, we notice that the range of the WD's mass decreases, while its average value 

increases with the cooling time. These results could provide a constraint when simulat- 
ing the SN la explosion, i.e. the WDs with a high C/O ratio usually have a lower central 
PsJ ■ density at ignition, while those having the highest central density at ignition generally 

K^ ' have a lower C/O ratio. The cooling time is mainly determined by the evolutionary age of 

Cn , secondaries, and the scatter of the cooling time decreases with the evolutionary age. Our 

^D ' results may indicate that WDs with a long cooling time have more uniform properties than 

!r^ , those with a short cooling time, which may be helpful to explain why SNe la in elliptical 

• ' galaxies have a more uniform maximum luminosity than those in spiral galaxies. 

o 
o 
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1 INTRODUCTION 



^ , As one of the most widely used distance indicators, type la supernovae (SNe la) show their importance 

j^ ■ in determining c osmological parameters, which resulted i n the discovery of the accelerating expansion 

of the universe dRiess et al.lll998l: iPerlmutter et al.|[l999l) . The result was exciting and suggested the 
presence of dark energy. At present, SNe la are proposed to be cosmological probes for testing the evo- 
lution of the dark energy equation of state with time and testing the evolutionary history of the universe 



JRiess et al.l 120071 : iKuznetsova et alj 20081: iHowell et al]|2009al) . They were even chosen to check the 
consistency of general relativity dZhao et alj|2010h . When SNe la are applied as a distance indicator, the 



Phillips relation is adopted, which is a linear relation between the absolute magnitude of SNe la at max- 
imum light and the magnitude drop of the B light curve during the first 15 days following the maximum 
(|Phillipslll993l) . This relation implies that the brightness of SNe la is mainly determined by one parame- 
ter. It is generally agreed that th e amount of ^ ^Ni formed during the supernova explosion dominates the 
maximum luminosity of SNe la ( Arnett 1982), but the origin of the different amount of ^^Ni for different 
SNe la is still unclear jPodsiadlowski et al...2008.) . Some numerical and synthetical resul ts showed that 
metallicity has an effect on the final amount of ^^Ni, and thus the maximum luminosity dTimmes et al.l 
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l2003UTravaglio et al.ll200l IPodsiadlowski et alll2006l: iBravo et al.ll2O10l) and there d o be some obser- 
vational evidence of the correlation between the properties of SNe la and metallicitv ([Branch & Berghl 
1 19931 lHamuvetal.lll996t IWang et al.|[l997l: ICappellaroet"aLlll997t IShanks et al.ll2002l) . However, the 
metallicitv see ms not to have the ab ilitv to interpret the scatter of the maximum lumi nosity of SNe la 
JTimmes et al.ll2003HGallagher et al.ll2008l: iHowell et ai]|2009bh . lN"omoto et alj ( fl999ll2003l) suggested 
that the ratio of carbon to oxygen (C/O) of a white dwarf at the moment of explosion is the dominant 
parameter for the Phillips relation. The higher the C/O, the larger the amount of nickel-56, and then the 
higher the maximum luminosity of SNe la. The C/O ratio is a function of the initial mass of the WD, 
which then is related to the progenitor system of the SNe la. By comparing theory and observations, the 
results of Meng et al. (2009) and Meng & Yang (2010a| upheld this suggestion. Lesaffre et al. (2006[) 
carried out a systematic study of the sensitivity of ignition conditions for H-rich Chandra single degener- 
ate exploders on various properties of the progenitors, and su ggested that the cent ral density of the WD 
at ignition may be the origin of the Phillips relation (see also IPodsiadlowski et al.ll20 08). These authors 
noticed that the more massive and/or the cooler the CO WD is when accretion begins, the higher the 
central density is at ignition. The central density is then also related to the progenitor system. When one 
simulates the explos ion of SNe la, the C/O and the central density are always set to be free parameters 
dRopke et al.ll2006al) . It is thus interesting to analyze whether there is a correlation between the C/O and 
the central density. In addition, analyzing how t he initial mass of the CO WD and its cooling time vary 
with the delay time is also an interesting task dGreggioluOlOl) . The purpose of this paper is to check 
these interesting problems. 

In section |2] we describe our model. We show the results in section [3] and give discussions and 
conclusions in sections ID 

2 MODEL AND PHYSICS INPUTS 

As suggested both iNomoto et"ai1 d 19991 l2003h and iLesaffre et"an (|2006|) . the progenitor model is the 
so-called single degenerate model (SD), i.e. the companion is probably a main seque nce or a shghtly 
evolved star (WD-hMS) or a red-giant (W D-hRG) or a helium star (WD + He star) (IWhelan & IbenI 



19731 : iNomoto. Thielemann & Yokoil 1984|). The SD model is widely accepted and studied by many 
authors dYungelson et alJll995l:lLi & van den Heuve]|l997l:lHachisu et al.lll999albl:lNomoto et alJll999l 
2003' 'Langer et al."2000': Han & Podsiadlowski'2004.'2006;'Chen & Li'2007'. '20091: lHanl2008l: lLu et al.l 



2009: Wang et al. 20 09a b: Wang, Li & Ha n 2010: Meng & Yang 2010b c). In this paper, we also ex- 
plore the SD model. IMeng & Yand d2010al) developed a compreh ensive progenitor mo del for SNe la. 
In the model, the mass-stripping effect by optic allv thick wind dHachisu et al.|[l996l) and the effect 



of a thermally unstab le disk were included dHachisu et al.ll20()8l : IXu & Lill2009l) . The prescription of 



or a tnermaiiy unstab le aisK were mciuaea QMacmsu et ai.iiz(j(j;a : lAu & LmzowA) . ine prescription or 
iHachisu et alj dl999al) for WDs ac creting hydrogen-rich m aterial from their companions was applied to 
calculate the WD mass growth. In lMeng & Yand d2010al) , both the WD + MS and WD + RG scenarios 
are considered, i.e. Roche lobe overflow (RLOF) begins at either the MS or the RG stage. After the 
RLOF, the WD accretes hydrogen-rich material from its donor and i ncreases its mass smoothly. When 
the mass of the WD reaches 1.378 Mq dNomoto, Thielemann & YokoL1984.) , the WD is assumed to 
explode as an SN la. They considered more than 1600 different WD close binary evolution and ob- 
tained a parameter space for SNe la which is summar ized in an orbital period - secondary mass plane 
(logP', Mj). Our exploration is based on the model of lMeng & Y ang (2010a). 

The C/O is a function of the initial WD mass of the progenitor system, i.e. a higher WD mass leads 
to a lower C/O. We then use the initial WD mass to represent the C/O. The central density is mainly 
determined by the initial WD mass and its cooling time, i.e. the more massive the WD and /or the longer 
the cooling time, the higher the central density at ignition. The cooling time of the CO WD is the time 
which elapses between its formation and the start of the accretion phase. We thus use the initial WD 
mass and its cooling time to represent the central density of a CO WD. The following is how we obtain 
the initial WD mass and the cooling time. 

To obtain the distributions of the initial WD mass and its co oling time, we carried ou t a series of 
Monte Carlo simulations via Hurley's rapid binary evolution code dHurlev et al.l200Clll2002h . The model 
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grids obtained by lMeng & YangI (1201 Oal) were incorporated into the code. We followed the evolution of 
10^ binaries and a circular orbit is assumed for all binaries. The metallicity is set to solar metallicity, 
i.e. Z = 0.02. The basic par ameters for the simulat ions are as follows: (1) a single star burst; (2) the 
initial mass function (IMF) of lMiller & Scald (119791) : (3) the mass-ratio distribution is constant; (4) the 
distribution of separations is constant in logo for wide binaries, where a is the orbital separation; (5) 
the common envelope (CE) ejection efficiency acE> which d enotes the fraction of th e released orbital 
energy used to eject the CE, is set to be either 1.0 or 3.0 (see lMeng & YangI (1201 Oah for details of the 
parameter input). 

There are three channels that produce WD + MS systems and one channel that produce a WD + 
RG system according to the situation of the primary in a primordial system at the onset of the first 
RLOF, i.e. the He star channel, the EAGB channel and the TPAGB channel (see lMeng & YandbOlol 
for details about the channels). The formed SD systems continue to evolve and the secondaries may fill 
their Roche lobes at a later stage, leading to the start of the RLOF. We assume that if the initial orbital 
period, P^^t^, and the initial secondary mass, Afj, of an SD system is located in the appropriate regions 
in the (log P', A/j) plane for SNe la at the onset of RLOF, an SN la is then produced. The cooling time 
is then the difference between the formation time of an SD system and the starting time of the RLOF. 



3 RESULTS 






^ d 



d 



^ 



CO 

d 



—1 1 1 1 1 1 1 1 1 r- 



—1 1 r- 



.''■■Jit' 



H 1- 



H 1 1 1- 



«CE=1° 



(1) 

H 1 1 1- 



~ T& 



«CE=3.0 




-:^^- 



(2) 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



9 

log(*oool) 



10 



Fig. 1 Distributions of initial WD masses and their cooling times for different acE- 



3.1 Distributions of the initial WD mass and their cooling time 



In Fig. [T] we show the distributions of the initial WD mass and its cooling times for different acE- We 
see from the figure that acE does not significantly affect the distributions. We do not find any correlation 
between the initial WD mass and its cooling time, which means that taking the C/O ratio and the central 
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Fig. 2 Distributions of the initial secondary mass and the cooUng times of CO WDs for 
different acE- 



density of a CO WD as free parameters is reasonable. However, the range of the WD mass decreases 
and the mean value of the WD mass incre ases with the cooling t ime, i.e. the range of the cooling time 
increases with the WD mass (see Fig. 9 in lMeng & Yandl2010al) . In addition, it is important to notice 
that there are some events with a very long cooling time, i.e. several Gyrs, althoug h the majority of CO 
WDs have a cooling time shorter than 1 Gyr. The models in lLesaffre et al.l (l2006h include cooling time 
ages of ~ 1 Gyr at most. For a long cooling time, the CO WD may become more degenerate before 
the start of the accretion phase, and some other pr ocess like C and O se paration or crystallization may 
occur, and dominate the properties of the CO WD (iFontaine et alj|200ll) . The simulation of the SNe la 
explosion under extremely degenerate conditions should then be a more important problem than that 
presented in literatures. 

We also investigate the distributions of the initial orbital period and the cooling time of a CO WD. 
Similar to the distribution of the initial WD mass and its cooling time, no correlation between the initial 
orbital period and the cooling time of a CO WD was discovered. 



3.2 Distributions of the initial secondary mass and the cooling time 

We present the distributions of the initial secondary mass and the cooling time for different acE values 
in Fig. |2] Although acE also does not significantly affect the distributions, these distributions differ 
remarkably from those of the initial WD mass and the cooling time. The cooling time is highly relevant 
to the secondary mass. It is clearly shown in Fig.|2]that there is an upper limit of the cooling time for a 
given secondary mass, and most of the events have a cooling time equal to the upper limit. The upper 
limit is derived from a constraint that the cooling time should be smaller than the evolutionary age of 
the secondary. However, there is still a scatter in the cooling time for a certain secondary mass, which 
originated from the different initial orbital period for an initial system with given WD and secondary 
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Fig. 3 Relation between the cooling time of CO WDs and the delay time of SNe la for 
different acE- 



masses. The shorter the initial orbital period, the earlier the onset of the RLOF between the WD and 
secondary, and then the shorter the cooling time. In addition, the scatter decreases with the decrease of 
secondary mass, and when the cooling time is larger than 2 Gyr, the scatter almost disappears. 



3.3 Relation between the cooling time of CO WDs and the delay time of SNe la 

It is interesting to analyze how the WD cooling time varies with the delay time of SNe la. In Fig. [3] the 
relation between the cooling time and the delay time is presented. The relation is very tight, and there 
is also an upper limit of the cooling time for a given delay time. This result is directly from the relation 
between the secondary mass and delay time. For most of the progenitor systems with a certain secondary, 
the delay time is the evolutionary lifetime of the secondary (Meng & Yang 2010a; Greggio 2010). So, 
using the secondary evolutionary time as the delay time of SNe l a is an excelle nt approximation when 
one applies an analytic method to study the delay time of SNe la (lGreggiol2010l) . There is also a scatter 
for the relation between the cooling time and the delay time, and the scatter decreases with the delay 
time. When the delay time is larger than 2 Gyr, the cooling time is almost equal to the delay time. 

4 DISCUSSION AND CONCLUSIONS 

In this paper, we do not find a correlation between the initial WD mass and its cooUng time. Since the 
C/O ratio is a function of the initial WD mass and the central density for a WD with giv en initial mass is 
mainly determined by the cooling time (iNomoto et alJl999ll2003HLesaffre et al.l2006l) . our results may 
imply that the C/O ratio and the central density at ignition are free parameters when simulating an SNe 
la explosion. However, our results still provide a constraint when simulating these types of explosions, 
i.e. the cooling time of WDs with an initial mass less than 1 AIq is generally shorter than 1 Gyr, but 
it may be as long as 15 Gyr for WDs with an initial mass larger than 1 AIq (see Fig. \T^. Because a 
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high initial WD mass me ans a lower C/O, and a massive WD and a long cooling time leads to a high 
central density at ignition dNomoto et al.lll999ll2003l:lLesaffre et al.ll2006) , the above result could imply 
that WDs with a high C/O usually have a lower central density at ignition, but those having the highest 
central density at ignition generally have a lower C/O. 

We also checked the effect of metallicity o n the distributions of the initial WD mass and the cooling 
time by Z ~ 0.001 (iMeng. Li & Yandl2010t) and no significant effect was found. Theory and obser- 
vations di d confirm t hat the effect of metallic ity cann ot explain the scatter in the maximum luminosity 
of SNe la (iTimmes et al. 2003: Gallagher et al. 2008; Howell et al. 2009b) . MetalUcity should, at most, 
be the s econdary parameter for the Phill ips relation, i.e. it is the origin for the scatter of the Phillips 
relation dPodsiadlowski et al.' 2006,7008). Furthermore, even regarding the secondary parameter, there 
is not a consensus. Mazzali & Podsiadlowski (2006) suggested that a variation of the relative content of 
(^^Fe+^^Ni) versus ^^Ni may be responsible for the observed scatter of the Phillips relation. However, 
iKasen. Ropke & WooslevI (l2009l) argued that the breaking of spherical symmetry is a critical factor in 
determining both the Phillips relation and the observed scatter around it. Then, the origin of both the 
Phillips relation and the scatter of the relation still should be investigated carefully. 

We found that the range of the initial WD mass decreases and the average WD mass increases with 
th e cooling time, which is similar to the relation between the initial WD mass and the delay time found 
by iMeng & Yand ( 1201 Oal) . In addition, the scatter of the cooling time also decreases with delay time. 
These results may indicate that the difference among the WDs with s short cooling time could be huge, 
but the properties of WDs with a long cooling time might be more uniform (see from Fig. [T]that all the 
WDs with a cooUng time longer than several Gyr have a mass larger than 1 .0 Mq). Since a long cooling 
time is equivalent to a long delay time (see Fig.[3]l, our results may imply that the properties of SNe la 
with long delay times might be more uniform than those with short delay times. It is widely known that 
there exists a scatter of the maximum luminosity of SNe la, and the scatter is affected by its environment. 
The most luminous SNe la always occur in spiral galaxies, but both spiral and elliptical galaxies are hosts 
for dimmer SNe la, which lead to a dimmer mean peak brightness in elliptical than i n spiral galaxies, i.e . 



the maximum luminosity of SNe la hosted in elliptical galaxies is more uni form (iHamuv et al. 1 9961) . 
These properties may be qualitatively interpreted by the model developed by iMeng & Yandj 20 1 Oal). at 



least if the C/O ratio is the origin of the Phillips relation. However, the quantitative study bv lRopke et al.l 
(l2006ah showed th at the C/O ha s only little impact on the amount of produced ^^Ni. It should be noticed 
that the results in Ropke et al.l (12006 ai) are model dependent, and are much different from those of 



iNomoto et al.l (Il999l l2003l) . Considering that the central density is mainly determined by the initial 



mass of the WD and its cooling time, the results in Figs. [T]and|3] i.e. the average values of the WD 
mass and the cooling time increase with the delay time, which means that the average of the central 
density increases with the del ay times of SNe la. A high central density at ignition leads to a larger 
amount of the produced ^^Ni (iRopke et al.ll2006al) . The effect of the ce ntral density on the amount of 
^^Ni seems then to be opposite with observations of lHamuv et al.l(ll996h . Furthermore, the variation of 
the amount of ^^Ni derived from the central density only amou nts to about 7%, whi ch can not be used 
to interpret the variation in the maximum luminosity of SNe la dRopke et al.ll2006al) . Perhaps, the C/O 
ratio, the c entral density and th e metallicity all contribute to the variation of the maximum luminosity 
of SNe la dRopke et al.l l2006b'). Then, which is the dominant parameter for the Phillips relation is still 
an open queation. In addition, as the cooling/delay ti me increases, the CO WDs become more and 
more degenerate, and even crystallization may occur dFontaine et al.l 120011) . What is the effect of the 
crystallization on the amount of ^^Ni should be an interesting problem. 
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